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Abstract 

In  this  work  we  present  a  model  of  the  near  field  plasma  plume  of  a  Pulsed  Plasma  Thruster  (PPT).  As  a 
working  example  we  consider  a  micro-PPT  developed  at  the  Air  Force  Research  Laboratory.  This  is  a 
miniaturized  design  of  the  axisymmetric  PPT  with  a  thrust  in  the  10  juN  range  that  utilizes  Teflon™™  as  a 
propellant.  The  plasma  plume  is  simulated  using  a  hybrid  fluid-PIC-DSMC  approach.  The  plasma  plume 
model  is  combined  with  Teflon™  ablation  and  plasma  generation  models  that  provide  boundary  conditions 
for  the  plume.  This  approach  provides  a  consistent  description  of  the  plasma  flow  from  the  surface  into  the 
near  plume.  The  magnetic  field  diffusion  into  the  plume  region  is  also  considered  and  plasma  acceleration 
by  the  electromagnetic  mechanism  is  studied.  Teflon™  ablation  and  plasma  generation  analyses  show  that 
the  Teflon™  surface  temperature  and  plasma  parameters  are  strongly  non-uniform  in  the  radial  direction. 
The  plasma  density  near  the  propellant  surface  peaks  at  about  1024  m"3  in  the  middle  of  the  propellant  face 
while  the  electron  temperature  peaks  at  about  4  eV  near  the  electrodes.  The  plume  simulation  shows  that  a 
dense  plasma  focus  is  developed  at  a  few  millimeters  from  the  thruster  exit  plane  at  the  axis.  This  plasma 
focus  exists  during  the  entire  pulse,  but  the  plasma  density  in  the  focus  decreases  from  about  2xl022  m‘3  at 
the  beginning  of  the  pulse  down  to  0.3xl022  in3  at  5  ps.  The  velocity  phase  is  centered  at  about  20  km/s  in 
the  axial  direction.  At  later  stages  of  the  pulse  there  are  two  ion  populations  with  positive  and  negative 
radial  velocity.  Electron  densities  predicted  by  the  plume  model  are  compared  with  near  field 
measurements  using  a  Herriot  Cell  technique  and  very  good  agreement  is  obtained. 
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Nomenclature 


Qj  -  Joule  heat 
Qr  -  radiation  heat 

Qf  -  heat  due  to  the  particle  convection 

Ne  -  electron  density 

Te  -  electron  temperature 

Vi  -  velocity  at  the  Knudsen  layer  edge 

T  -  ablation  rate,  [kg/m2s] 

Ni,  N2  -  densities 

m  -  heavy  particle  mass 

Tb  T2  -  temperatures 

Cs  -  sound  speed 

j  -  current  density 

B  -  magnetic  field 

E  -  electric  field 

me  -  electron  mass 

Ve  -  electron  velocity 

vei  -  electron-ion  collision  frequency 

a  -  plasma  (Spitzer)  conductivity 

p  -  permittivity 


cor  -  Hall  parameter 


1.  Introduction 


The  pulsed  plasma  thruster  (PPT)  was  among  the  first  of  various  electrical  propulsion  concepts 
accepted  for  space  flight  mainly  due  to  its  simplicity  and  hence  high  reliability1.  However,  the 
PPT  has  an  efficiency  that  is  generally  low2  at  about  10%  leaving  open  the  opportunity  for 
considerable  improvement3.  Currently,  PPT’s  are  considered  as  an  attractive  propulsion  option 

for  stationkeeping  and  drag  makeup  purposes  on  mass  and  power-limited  satellites4,5. 
Guaranteeing  successful  operation  of  spacecraft  using  a  PPT  requires  a  complete  assessment  of 
the  spacecraft  integration  effects.  The  PPT  plume  contains  various  ion  and  neutral  species  due  to 
propellant  decomposition  and  possible  electrode  erosion.  Some  attempts  of  PPT  plume  modeling 
using  particle  simulations  were  performed  recently6,7,8.  In  Refs.  7  and  8  we  have  considered  the 
plume  flowfield  exhaust  from  the  recently  developed  electrothermal  PPT  (so-called  PPT-4)  and 
therefore  electromagnetic  effects  in  the  plume  were  neglected.  Different  variations  of 
electromagnetic  PPTs  are  also  candidates  for  various  missions9.  Recently,  a  micro-PPT  has  being 
designed  at  the  Air  Force  Research  Laboratory  (AFRL)  for  delivery  of  very  small  impulse  bit10. 
This  is  a  simplified  miniaturized  version  of  a  conventional  PPT  designed  to  provide  attitude 
control  and  stationkeeping  for  microsatellites.  We  will  use  the  AFRL  micro-PPT  as  a  working 
example  for  several  reasons.  Firstly,  electromagnetic  (jxB)  acceleration  is  the  primary 
mechanism  in  this  thruster;  and  secondly,  there  is  no  internal  flow  in  this  device  and  therefore  the 
near-field  plasma  plume  is  an  essential  part  of  the  thrust  generation  process.  Therefore  careful 
modeling  of  the  acceleration  is  needed  to  understand  the  characteristics  of  the  device  as  a  whole 
m  addition  to  being  a  pre-cursor  to  accurate  estimation  of  contamination  issues.  Since  in  this 
device  there  is  no  separation  between  the  main  plasma  acceleration  region  and  the  plume 
expansion,  both  regions  must  be  simulated  in  one  model.  Because  the  plasma  acceleration  is 
external,  the  plasma  is  sufficiently  rarefied  so  that  an  Magneto  hydrodynamic  (MHD)  code  such 
as  MACH2  (Ref.  1 1)  cannot  be  used  in  end-to-end  simulation. 
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An  accurate  model  of  the  PPT  plume  relies  on  the  boundary  and  initial  conditions.  These 
conditions  can  be  formulated  by  consideration  of  the  Teflon™  ablation  process.  The  Teflon™ 
ablation  computation  is  based  on  a  recently  developed  kinetic  ablation  model12,13.  In  this  model, 
the  detailed  physics  of  the  Teflon™  evaporation  is  studied  by  consideration  of  the  distribution 
function  of  the  particles  in  the  kinetic  layer  adjacent  to  the  surface. 

Another  important  effect  related  to  the  plasma  plume  exhaust  from  an  electromagnetic  PPT  is  the 
•magnetic  field  diffusion  into  the  near  plume.  Previously,  we  have  modeled  the  effect  of  the 
magnetic  field  on  the  near-field  plume  for  Hall  thrusters14  under  steady  state  conditions.  It  was 
found  that  the  magnitude  of  the  magnetic  field  at  the  thruster  exit  has  an  important  effect  on  the 
plasma  potential  distribution  in  the  plume.  In  the  present  research,  it  is  proposed  to  include  the 
electromagnetic  effects  on  the  near  field  plume  of  unsteady  plasma  flow.  The  computational 
domain  is  shown  in  Fig.  1 .  The  model  is  based  on  a  hybrid  approach  involving  a  direct  simulation 
Monte-Carlo  (DSMC)  description  of  neutrals,  a  Particle-In-Cell  (PIC)  model  for  ions,  and  a  fluid 
description  of  the  electrons.  In  these  methods,  the  potential  distribution  is  usually  calculated  by 
reducing  the  electron  momentum  equation  to  the  Boltzmann  relation  in  the  absence  of  a  magnetic 
field.  In  the  plasma  plume  domain  where  the  magnetic  field  exists,  i.e.  the  near  field  plume 
region,  it  is  necessary  to  include  the  magnetic  field  effects  in  the  electron  momentum  equation. 

2.  The  model  of  the  plasma  layer 

C 

The  model  presented  here  describes  the  plasma  layer  near  the  Teflon™  surface  as  shown  in  Fig. 

2.  The  model  of  the  plasma  layer  includes  Joule  heating  of  the  plasma,  heat  transfer  to  the 
Teflon™,  and  Teflon™  ablation.  Mechanisms  of  energy  transfer  from  the  plasma  column  to  the 
wall  of  the  Teflon™  include  heat  transfer  by  particle  convection  and  by  radiation.  The  Teflon™ 
ablation  computation  is  based  on  a  recently  developed  kinetic  ablation  model12  (see  next  section). 
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It  is  assumed  that  within  the  plasma  layer  all  parameters  vary  in  the  radial  direction,  r.  The  energy 
balance  equation  can  be  written  in  the  form: 


2Ne  dTe/dt  -  Qj  -  Qr  -  Qf- 


(1) 


This  equation  depends  on  the  coordinate  along  the  propellant  face.  For  known  plasma  density  and 
temperature  the  heat  flux  to  the  surface  is  calculated.  The  Teflon™  surface  temperature  is 
calculated  from  the  heat  transfer  equation  with  boundary  conditions  that  take  into  account 
vaporization  heat  and  conductivity.  The  solution  of  this  equation  is  considered  for  two  limiting 
cases  of  substantial  and  small  ablation  rate  very  similar  to  that  described  in  Ref.  8.  The  density  at 
the  Teflon™  surface  is  calculated  using  the  equilibrium  pressure  for  Teflon™.  The  plasma 
density  in  the  layer  is  determined  in  the  framework  of  the  kinetic  ablation  model  (see  next 
section).  For  known  pressure  and  electron  temperature  one  can  calculate  the  chemical  plasma 
composition  assuming  LTE8’15’16.  The  Saha  equations  are  supplemented  by  the  conservation  of 
nuclei  and  quasi-neutrality. 

3.  Ablation  model 

The  Teflon™  ablation  is  modeled  in  the  framework  of  the  approximation13  based  on  a  kinetic 
model  of  the  material  evaporation  into  discharge  plasmas12.  The  model  couples  two  different 
layers  between  the  surface  and  the  plasma  bulk  as  shown  in  Fig.  2:  (1)  a  kinetic  non-equilibrium 
layer  adjusted  to  the  surface  with  a  thickness  of  about  one  mean  free  path;  and  (2)  a  collision- 
dominated  layer  with  thermal  and  ionization  non-equilibrium.  The  velocity  at  the  edge  of  the 
kinetic  layer,  V1;  can  be  determined  from  the  coupling  solution  of  the  hydrodynamic  layer  and  the 
quasi-neutral  plasma.  For  known  velocity  and  density  at  this  interface,  it  is  possible  to  calculate 
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the  ablation  rate.  In  the  hydrodynamic  layer  the  relation  between  the  velocities,  temperatures  and 
densities  at  the  boundaries  1  and  2  as  well  as  the  ablation  rate  are  formulated  according  to  Ref.  13 
in  the  form: 


r  =  mV,N1  =N,[(2kT1/m)-(T2N2/2T1  -^/(N.-N,2/^)]0'5 . (2) 

The  system  of  equations  is  closed  if  the  equilibrium  vapor  pressure  can  be  specified  that 
determines  parameters  (N0  and  T0)  at  the  Teflon™  surface.  The  solution  of  the  Knudsen  layer 
problem  relates  parameters  at  the  boundary  1  to  the  parameters  at  the  boundary  0  (Ref.  12).  The 
full  self-consistent  solution  of  this  problem  can  be  obtained  when  the  ablation  is  coupled  with  the 
plasma  plume  expansion.  In  the  present  work  in  order  to  simplify  the  problem,  we  will  assume 
that  the  plasma  accelerates  up  to  the  sound  speed  near  the  boundary  2.  This  assumption  can  be 
justified  by  the  fact  that  due  to  significant  electrodynamic  acceleration  in  this  type  of  PPT,  the 
plasma  density  will  quickly  decrease,  therefore  providing  solution  of  the  ablation  problem  close 
to  that  for  ablation  into  vacuum.  In  this  case  the  plasma  density  at  the  edge  of  the  kinetic  layer 
will  be  equal  to  0.34-No  and  the  temperature  is  0.7-To.  The  flux  returned  to  the  surface  is  equal  to 
16%  of  the  ablated  flux  (Ref.  12). 

4.  Plasma  plume  electrodynamics 

The  general  approach  for  the  plume  model  is  based  on  a  hybrid  fluid-particle  approach  that  was 
used  previously  (Ref.  7).  In  this  model,  the  neutrals  and  ions  are  modeled  as  particles  while 
electrons  are  treated  as  a  fluid.  Elastic  (momentum  transfer)  and  non-elastic  (charge  exchange) 
collisions  are  included  in  the  model.  The  grids  employed  in  this  computation  are  also  similar  to 
those  used  previously  (Ref.7).  The  particle  collisions  are  calculated  using  the  DSMC  method.17 
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Momentum  exchange  cross  sections  use  the  model  of  Dalgamo  et  aln,  while  charge  exchange 
processes  use  the  cross  sections  proposed  by  Sakabe  and  Izawa.19  Acceleration  of  the  charged 
particles  is  computed  using  the  PIC  method.20  The  plasma  velocity  distribution  depends  upon  the 
magnetic  field  distribution  and  ion  dynamics  is  calculated  as  follows: 

dV/dr  =  -Cs2Vln(N)  +  jxB/mN . . (3) 

The  electron  dynamics  is  very  important  in  the  plasma  plume.  Previously  our  model  was  based  on 
the  assumption  that  electrons  rapidly  reach  the  equilibrium  distribution  and  in  the  absence  of  the 
magnetic  field  can  be  described  according  to  the  Boltzmann  distribution.  While  this  was  a 
satisfactory  assumption  in  the  case  of  an  electrothermal  thruster  plume  this  is  not  suitable  for  the 
near  field  of  an  electromagnetic  thruster.  In  the  presence  of  a  strong  magnetic  field,  the  electron 
density  distribution  deviates  from  that  according  to  Boltzmann21.  In  the  case  of  a  magnetic  field 
the  electron  momentum  equation  reads  (neglecting  electron  inertia): 

0  =  -e2Ne(E+VexB)  -  eVPe  -  veimj . (4) 

We  have  assumed  quasi-neutrality  therefore  Ne  =  N;  =  N.  The  electric  and  magnetic  field 
distributions  in  the  plume  are  calculated  from  the  set  of  Maxwell  equations.  We  further  assume 
that  the  magnetic  field  has  only  an  azimuthal  component  and  also  neglect  the  displacement 
current.  The  combination  of  the  Maxwell  equations  and  electron  momentum  conservation  gives 
the  following  equation  for  the  magnetic  field: 

dBfdt  =  l/(ap) V2B  -  Vx(jxB/(eN))  +  Vx(VxB) . .  (5) 
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A  scaling  analysis  shows  that  the  various  terms  on  the  right  hand  side  of  Eq.  5  may  have 
importance  in  different  regions  of  the  plasma  plume  and  therefore  a  general  end-to-end  plasma 
plume  analysis  requires  keeping  all  terms  in  the  equation.  In  the  case  of  the  near  plume  of  the 
micro-PPT  with  a  characteristic  scale  length  of  about  1  cm  the  magnetic  Reynolds  number 
Rem«l  and  therefore  the  last  term  can  be  neglected.  Taking  this  into  account  in  dimensionless 
form,  Eq.  5  can  be  written  as: 

Rem9B/5t  =  V2B  -  (cox)-  {Vx(VxBxB)} . (6) 

where  (got)  is  the  Hall  parameter  that  measures  the  Hall  effect.  Therefore,  depending  on  the 
plasma  density,  the  Hall  effect  may  be  important  for  the  magnetic  field  evolution.  One  of  the  first 
calculations  of  the  plasma  flow  with  Hall  effect  were  performed  by  Brushlinski  and  Morozov  (see 
Ref.  22  and  references  therein).  They  considered  isothermal  flow.  The  plasma  density  becomes 
high  at  the  cathode  and  lower  at  the  anode.  The  Hall  effect  has  a  particularly  noticeable  influence 
on  the  magnetic  field  distribution.  The  field  near  the  anode  increases  and  near  the  cathode 
decreases.  As  a  result  the  current  is  deflected  to  the  side  and  grazes  the  anode. 

Our  estimations  show  that  the  Hall  parameter  cor«l  if  the  plasma  density  near  the  Teflon™ 
surface  N>1023  m'3.  This  case  is  realized  in  the  micro-PPT  (see  the  next  section)  so  the  Hall  effect 
is  expected  to  be  small  for  this  particular  case.  Therefore  all  results  presented  below  are 
calculated  without  considering  Hall  effect.  Having  the  magnetic  field  distribution  one  can 
calculate  the  current  density  distribution  from  Ampere’s  law: 

Jlj  =  VxB . . . . . . . ; . (7) 


8 


The  magnetic  field  and  current  distributions  calculated  from  this  model  are  used  in  PIC  to 
evaluate  the  ion  dynamics  according  to  Eq.  3 . 

5.  Boundary  conditions 

The  boundary  conditions  for  the  magnetic  field  calculations  are  shown  in  Fig.  1 .  We  assume  that 
the  current  is  uniform  on  both  electrodes  that  allows  us  to  estimate  the  current  density  on  the 
cathode  jc  and  on  the  anode  ja.  The  magnetic  field  is  assumed  to  vary  as  1/r  on  the  upstream 
boundary.  At  the  lateral  boundary  we  assume  that  the  normal  current  jn=0.  The  downstream 
boundary  is  considered  to  be  far  enough  away  that  B=0  can  be  assumed.  Along  the  centerline  the 
magnetic  field  is  zero. 

The  boundary  conditions  for  the  plume  are  generated  through  solution  of  the  Teflon™  ablation 
problem  as  will  be  presented  in  the  Results  section.  These  involve  time  and  radial  dependent 
variations  of  the  plasma  (including  Carbon  and  Flourine  ions  and  neutrals)  density  and  electron 
temperature. 

The  results  are  presented  for  a  3.6  mm  (0.141”)  diameter  micro-PPT,  which  has  a  0.9  mm 
diameter  central  electrode,  3.1  mm  propellant  diameter  and  0.24  mm  anode  wall  (Ref.  10).  In 
these  simulations,  the  experimental  current  waveform  is  used,  that  is  described  in  a  first 
approximation  as  an  underdamped  LRC  circuit  current: 

I(t)  =  Ip-sin(at)exp(-pt) . (8) 
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where  Ip  = 


L  is  the  effective  inductance  in  the  circuit,  C  is 


the  capacitance,  R  is  the  total  circuit  resistance,  and  E  is  the  pulse  energy.  Results  presented 
below  correspond  to  E=2.25  J  and  C=0.5  pF.  The  best  fit  with  the  experimental  waveform 
(frequency)  corresponds  to  a=4.7-107  rad/s  and  circuit  inductance  L=90  nH. 


6.  Results 


The  plasma  density  and  electron  temperature  distribution  are  also  shown  in  Figs.  3a-3b.  The 
plasma  density  peaks  at  about  1024  m'3  midway  between  the  electrodes.  The  electron  temperature 
is  strongly  non-uniform  radially  with  peaks  near  the  electrodes  of  about  4.5  eV  as  shown  in  Fig. 
3b.  The  reason  for  higher  electron  temperature  near  the  electrodes  is  due  to  current  spreading  in 
the  space  between  the  electrodes  and  current  focusing  near  the  electrodes  (see  below  results  on 
current  distribution). 

The  spatial  and  temporal  variation  of  the  Teflon™  surface  temperature  is  shown  in  Fig.  3c.  The 
Teflon™  temperature  sharply  increases  during  the  first  2  ps  of  the  pulse  and  peaks  at  about  960 
K.  One  can  see  that  the  temperature  is  generally  non-uniform  in  the  radial  direction  and  has  a 
minimum  at  radial  distances  of  1.1 -1.3  mm.  Since  the  Teflon™  ablation  is  approximately 
exponentially  proportional  to  the  surface  temperature,  the  model  predicts  a  lower  rate  of  ablation 
in  the  areas  where  the  surface  temperature  has  a  minimum.  Taking  this  into  account,  the  effect  of 
the  temperature  distribution  may  be  related  to  the  preferential  charring  of  the  Teflon™  surface 
observed  experimentally  [Ref.  10].  As  was  mentioned  earlier,  the  ablation  rate  is  also  non- 
uniform  radially.  This  effect  is  shown  in  Fig.  3d.  One  can  see  that  the  ablation  rate  peaks  near  the 
electrodes  at  about  60  kg/m2s,  while  in  the  middle  of  the  propellant  face  it'  is  about  30-40  kg/m2s. 
The  calculated  total  ablated  mass  per  pulse  is  about  1 .4  pg. 
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A  region  of  magnetic  field  diffusion  in  the  near  field  outside  the  micro-PPT  is  shown  in  Fig.  4a. 
The  magnetic  field  drops  by  an  order  of  magnitude  at  about  1 .5  mm  that  is  equal  to  the  thruster 
radius.  This  is  also  the  region  where  most  of  the  current  is  concentrated  as  shown  in  Fig.  4b.  One 
can  see  that  the  current  density  is  high  near  the  central  electrode  and  near  the  outer  electrode.  This 
is  a  reason  for  the  increasing  Teflon™. surface  temperature  and  electron  temperature  in  these 
regions.  According  to  the  model  presented  in  Sec.  4  the  electromagnetic  acceleration  of  the 
plasma  is  also  expected  to  occur  in  this  region. 

Figure  5  shows  evolution  of  the  Carbon  ion  (C+)  component  of  the  plasma  plume  during  the  main 
part  of  the  pulse.  One  can  see  that  a  dense  plasma  focus  is  developed  at  a  few  millimeters  from 
the  thruster  exit  plane.  This  plasma  focus  exists  during  the  entire  pulse  as  shown  in  Fig.  5,  but  the 
plasma  density  in  the  focus  decreases  from  about  2xl022  m'3  at  the  beginning  of  the  pulse  down  to 
0.3xl022  m'3  at  5  ps.  At  the  beginning  (first  2  ps)  the  C+  density  mainly  develops  a  gradient  in 
the  radial  direction  that  is  a  result  of  high  directed  velocity  in  the  axial  direction.  Later,  during  the 
pulse,  the  axial  density  gradient  becomes  comparable  to  the  radial  one. 

The  Flourine  ions  (F+),  due  to  their  larger  mass,  have  different  dynamics  as  shown  in  Fig.6.  They 
have  smaller  acceleration  in  the  axial  direction  even  at  the  beginning  of  the  pulse  and  therefore 
both  axial  and  radial  density  gradients  are  developed.  The  F+  density  in  the  plume  and  in  the 
plasma  focus  is  larger  than  that  of  C+,  because  originally  Teflon™  has  composition  C2F4  with  F 
density  twice  larger  than  that  of  C.  Additionally  F  ions  experience  less  acceleration  in  the  plume 
because  of  their  higher  mass  that  also  contributes  to  their  relative  density  increase. 
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The  micro-PPT  is  essentially  an  electromagnetic  accelerator  as  shown  in  the  velocity  phase  plots 
'  (Figs.  7,8).  The  phase  plot  of  the  Carbon  ions  at  1  ps  is  centered  at  20  km/s  in  the  axial  direction. 
Ions  experience  also  radial  expansion  in  both  directions  due  to  the  magnetic  field  structure  and 
the  temperature  expansion.  The  radial  velocity  in  the  negative  direction  is  related  to  the  focus 
formation  along  the  axis,  as  shown  in  Figs.  5,  6.  The  Flourine  ions  generally  have  both  smaller 
axial  and  radial  velocities  due  to  their  higher  mass.  At  a  later  stage  of  the  pulse  (see  Fig.  8) 
clearly  there  are  two  ion  populations  with  positive  and  negative  radial  velocities.  This  is  due  to 
the  annular  plasma  injection  corresponding  to  the  thruster  geometry  (see  Figs.  1,2). 

During  the  entire  pulse  there  is  a  population  of  ions  having  a  negative  axial  velocity  with  the 
magnitude  up  to  about  10  km/s  (see  Figs.  7,8).  This  population  creates  the  backflow 
contamination  that  occurs  mainly  onto  the  thruster  itself.  The  Carbon  ions  have  a  larger  negative 
velocity  due  to  their  higher  mobility  that  results  in  their  domination  in  the  backflux.  This  backflux 
may  be  mainly  responsible  for  charring  phenomena  observed  in  this  thruster  (see  Ref.  10). 

7.  Comparison  with  experiment 

In  this  section  we  present  measured  and  predicted  electron  density  distributions  in  the  near  field 
plume  for  one  micro-PPT  design.  These  data  will  be  compared  in  order  to  assess  our  plume  and 
device  model. 

An  experimental  basis  for  comparison  is  provided  using  a  Herriott  Cell  interferometer.  Electron 
density  measurements  are  taken  on  a  6.35  mm  (1/4”)  diameter  micro-PPT  at  AFRL.  The 
interferometer  uses  a  single  laser  wavelength  and  quadrature  heterodyne  technique  described  by 
Spanjers  et  al23  Addition  of  a  Herriott  Cell  acts  to  confine  a  large  number  of  laser  passes  into  an 
area  suitable  for  maximum  exposure  to  the  MicroPPT  plume.  This  is  achieved  by  focusing  the 
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laser  between  the  two  concave  mirrors  of  the  cell.  The  technique  is  used  to  increase  signal-to- 
noise  ratio  for  diffuse  plasmas  by  increasing  laser  exposure  to  the  plasma  over  a  characteristic 
path  length.24  Thirteen  laser  reflections  in  the  Herriott  Cell  were  focused  to  two  points,  separated 
by  3  mm.  For  data  shown  here,  these  points  formed  a  plane  parallel  to  the  fuel  face  and  5  mm 
distant.  More  details  about  the  Herriott  Cell  technique  can  be  found  elsewhere.24 

Figure  9  shows  the  experimental  data  co-plotted  with  model  predictions.  Plasma  density  peaks  at 
about  3x1  O' 6  cm*3  and  decreases  by  several  orders  of  magnitude  towards  the  pulse  end.  The 
experimental  data  was  taken  at  a  discharge  energy  of  6.6  J  from  a  0.417  pF  capacitor. 
Experimental  waveforms  of  the  current  were  obtained  using  a  self-integrating  Rogowski  coil. 
Peak  density  reaches  23±6T015  cm'3  with  uncertainty  due  to  shot-to-shot  variations  in  thruster 
firing.  One  can  see  that  the  model  correctly  predicts  both  the  plasma  density  level  and  temporal 
behavior  during  the  entire  pulse. 
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7.  Concluding  remarks 


In  this  paper,  a  self-consistent  description  of  an  electromagnetic  pulsed  plasma  thruster  from 
plasma  generation  into  the  near  plume  was  presented.  A  micro-PPT  developed  at  AFRL  was 
considered  as  a  working  example.  In  this  device,  no  separation  exists  between  the  main  plasma 
acceleration  region,  which  usually  occurs  in  an  internal  flow,  and  the  external  plasma  plume  field. 
Therefore,  a  single  end-to-end  model  is  necessary  for  accurate  simulations.  A  kinetic  Teflon™ 
ablation  model  was  incorporated  in  order  to  provide  the  boundary  conditions  for  the  plasma 
plume. 

The  phenomena  in  the  plasma  plume  related  to  the  electromagnetic  effects  were  studied.  The 
plume  simulation  showed  that  a  dense  plasma  focus  developed  at  a  few  millimeters  from  the 
thruster  exit  plane  at  the  axis.  This  plasma  focus  exists  during  the  entire  pulse,  but  the  plasma 
density  in  the  focus  decreases  from  about  2xl022  m'3  at  the  beginning  of  the  pulse  down  0.3xl022 
m’3  at  5  j_ls.  The  velocity  phase  is  centered  at  about  30  km/s  in  the  axial  direction  demonstrating 
that  the  micro-PPT  is  essentially  an  electromagnetic  accelerator.  At  a  later  stage  of  the  pulse  there 
are  two  ion  populations  with  positive  and  negative  radial  velocity.  It  is  predicted  that  there  is  a 
population  of  ions  having  a  negative  axial  velocity  magnitude  up  to  about  10  km/s.  This 
population  creates  the  backflow  contamination  that  flows  mainly  onto  the  thruster  itself.  The 
Carbon  ions  have  a  larger  negative  velocity  due  to  their  higher  mobility  that  results  in  their 
domination  in  the  backflux.  It  is  believed  that  this  backflux  is  responsible  for  the  charring 
phenomena  observed  in  this  thruster. 

Predicted  electron  density  in  the  near-field  plume  was  directly  compared  with  experimental  data 
and  very  good  agreement  was  obtained. 
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Figure  Captions 


Figure  1.  Schematic  diagram  of  micro-PPT  plume  and  boundary  conditions 


Figure  2.  Schematic  of  the  near  Teflon™  plasma  layer 


Figure  3.  Spatio-temporal  distribution:  (a)  Teflon™  surface  temperature,  (b)  plasma  density,  (c) 
electron  temperature,  and  (d)  ablation  rate. 


Figure  4.  (a)  Magnetic  field  distribution  and  (b)  current  line  vectors  in  the  near  field  of  the  micro- 
PPT. 


Figure  5.  Evolution  of  the  Carbon  ion  density  during  the  pulse 


Figure  6.  Evolution  of  the  Flourine  ion  density  during  the  pulse 


Figure  7.  Ion  velocity  phase.  Early  stage  of  the  pulse 


Figure  8.  Ion  velocity  phase.  Late  stage  of  the  pulse 


Figure  9.  Comparison  of  predicted  and  measured  electron  density  time  variation  at  5  mm  from 
the  propellant  face  at  the  axis  in  the  case  of  the  6.35  mm  diameter  micro-PPT  firing  at  6.6  J. 
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Figure 2.  Keidar  et  al.  “Electromagnetic  effects...” 
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Figure  4.  Keidar  et  al.  “Electromagnetic  effects 
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Figure  6.  Keidar  et  al.  “Electromagnetic  effects. . 
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Figure  7.  Keidar  et  al.  “Electromagnetic  effects 
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Figure  8.  Keidar  et  al.  “Electromagnetic  effects. . . 
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